Correlation energy of two electrons in the high-density limit 
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We consider the high-density-limit correlation energy E c in D > 2 dimensions for the S ground 
states of three two-electron systems: helium (in which the electrons move in a Coulombic field), 
spherium (in which they move on the surface of a sphere), and hookium (in which they move in a 
quadratic potential). We find that the E c values are strikingly similar, depending strongly on D 
but only weakly on the external potential. We conjecture that, for large D, the limiting correlation 
energy in any confining external potential, where S = 1/(D — 1). 
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I. INTRODUCTION 



The concept of electron correlation energy (E c ) is an 
old and important one, first introduced by Wigner [l| 
and later defined by Lowdin 0] as the error 



E c — E — E 



HF 



(1) 



of the Hartree-Fock (HF) model. Understanding and cal- 
culating E c is one of the most important and difficult 
problems in quantum chemistry and molecular physics. 

The observation that HF theory is useful for the pre- 
diction of molecular structure |3| suggests that E c of- 
ten depends only weakly on the external potential. To 
explore this, we have studied the correlation energy 
E c (D,m, Z) of two opposite-spin electrons, confined to 
a D-dimensional space and moving in an external poten- 
tial Z m+2 V(r) where V(r) oc r m . 

We consider three such systems. In D-helium, the elec- 
trons move in the Coulomb potential V(r) = — 1/r. In D- 
spherium jij, they move in a constant potential V(r) = r° 
on a £>-sphere of radius l/Z. (In 2-sphcrium, for exam- 
ple, this is the surface of a three-dimensional ball.) In 
_D-hookium (also known as Hooke's atom or harmonium) 
they move in the harmonic potential V(r) = r 2 /2. 

After the length scaling r <— Zr, the Hamiltonians of 
the three systems reduce to the form [|| 0] 



Many workers have investigated the energies 

EJD,m)= lim E c = E 2 (D,m) - Ef F (D,m) (5) 

Z— too 

that arise in the high-density limit. Studies of the helium- 
like ions @, Q, for example, showed that E c (3, —1) ps 47 
mh (millihartrees) and several groups have noted that 
the limits for 3-spherium \E C (3,0) « 48 mh] (T^j and 
3-hookium [£^(3, 2) w 50 mh] (Tm3] are similar. 

Two-dimensional systems have also been studied and, 
although the limiting energies for 2-hclium \E C (2, — 1) > 
212 mh] [3, 2-spherium [^ c (2,0) « 227 mh] [H El] 
and 2-hookium [E c (2, 2) > 162 mh] [l7| are several times 
greater than their D = 3 analogs, they appear similar to 
one another. Conversely, those for 4-helium [EJA, — 1) > 
18 mh] and higher heliums are much smaller (14| . 

Such results suggest that the limiting correlation en- 
ergies arc similar, not only for D = 3 as previously re- 
ported, but also for other D. It leads to the idea that 
the correlation energy of two electrons in the high-density 
limit depends strongly on the dimensionality of the space 
in which they move, but weakly on the external potential. 
To explore this, we have calculated the limiting correla- 
tion energies of the 1 S ground states of helium, spherium 
and hookium for D = 2,3, ... ,8. We use atomic units 
throughout. 



X 
S3 



H 



V? 



^i + V(ri) + V{r 2 ) + -±- (2) 



where r\2 = \r\ — r%\. Following Hylleraas [8j, perturba- 
tion theory can then be used to expand both the exact 
and HF energies as series in l/Z, yielding 

E = E {D,m)Z 2 + E 1 (D,m)Z + E 2 (D,m) + OiZ- 1 ) 

(3) 

Ehf = E (D, m)Z 2 + Ei(D, m)Z + Ef F {D, m) + O^ 1 ) 

(4) 
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II. HELIUM 



The one-electron Hamiltonian in D-helium is 

d 2 D — 1 d 

dr 2 



o 



r dr 

and the zeroth-ordcr wave function is 



*o(r-i,r 2 ) = 



(D-1) D T(D) 
The Eq and E\ values are 

4 



exp 



2n + 2r 2 
D-l 



E (D,-1) 
Ei(D, — 1) 



(D-l) 2 
4 



T(D+l)T(S±l) 



(D-i) 2 r(D + i)r(f) 



(6) 

(7) 

(8) 
(9) 



2 



where T is the Gamma function [181 ]. 

Ei values were computed using the Hylleraas method 
[H . We adopted the length and energy scaling of Herrick 
and Stillinger [l9| and used the Hylleraas basis functions 



where, following previous work [28|, |2S| . we use 8 
1/p-l). 



III. SPHERIUM 



s n t l u m exp(-s/2) 



(10) 



where s = r\ + ?*2 , t = r\ — r 2 and u = 7*12 • The second- 
order energy, which minimizes the Hylleraas functional, 
is then given by 



E 2 (D,-1) 



ib T A- x b 

2 



(11) 



where 



D - 1 



-2E (12) 



b w — 2EiSquj t: — Vq w 



(13) 



with uj = (n,2l,m). In (12]) and (T3|), M, L, S and V 
are the kinetic, electron- nucleus, overlap and repulsion 
matrices, respectively. Details can be found elsewhere 
fl9l |20T ] . The Hylleraas basis was progressively enlarged 
by increasing the maximum values of n, I and m until 
the most difficult case (D — 2) converged to 6 digits. 

Although the E 2 value for 3-hclium has been studied in 
great detail (as in, for example, the work of Morgan and 
co-workers ({J UJ ) , the only other helium whose E2 value 
has been reported (l9| (by exploiting interdimensional 
degeneracies [22|, [23|) is 5-helium. 

Although Loeser and Herschbach have investigated the 
dimensional dependence of the HF energy of helium 0, 
[H, Ef F (D, -1) has been reported [H only for D = 3. 
All values can be found using the generalization 



E™(D,-1) = - f 
Jo 

W{r) = 2 



W{rf 



„D-1 



dr 



(14) 



'- -*oM 2 

[J{x)-E 1 ]^ {xfx D - l dx (15) 
of the Byers-Brown-Hirschfelder equations (26j , where 



J(r) 



max(r, x) 



-F 



3-D 1 D 



2 ' 2' 2 



, a 



x^dx 



min(a;,r) 
max(x,r) 

For odd D, this yields simple expressions such as 



(16) 



:v and F is the hypergeometric function [1! 



£ 2 HF (3,-1) 
£ 2 HF (5,-1) 



9 1 3 13 
— In (17) 

32 4 432 v ' 

903 , 3 35 213 

In - - , _ , — (18) 



1024 4 124 416 
Eqs (|TT|) and (fl4|) yield the large- D expansions [1 



E 2 (D,-l)~--5 



5 -2 31 ^3 



384 



E™(D,-1) 



-l S 2 + —8U 
2 32 

v ' ; 8 384 



(19) 
(20) 
(21) 



The zeroth-order Hamiltonian of D-spherium is 

H =~-(D-l)cat9i (22) 



(where 9 is the inter-electronic angle) and the associated 
eigenfunctions and eigenvalues are, respectively, 



tf B (0) =ATC n 2 (cos 6) 
e„ = n(n + D - 1) 

where C n 2 is a Gegenbauer polynomial [l8| and 



<2 D - 3 (2n + D - i)r(£=!) T(n+1 



7rr(?i + D- 1) 



Using the partial- wave expansion of , one finds 



12 



£jl \ _ (n + l) D - 2 



(n+|) 



(23) 
(24) 



(25) 



(26) 



where (a) n is a Pochhammer symbol [la ] and, therefore, 

r(£>-i)r(^) 



r(D-i)r(f) 

The second-order energy is given by 

(^ol^l^n) 2 



(27) 



E 2 (D,0) = ^2 



£0 — £n 



r(D)r(V 



-in 



E 



(n + 1) 



D-2 



r(f) 2 ^(n + i)^ 



1 



1 

n n + D — I 
(28) 



which reduces to gen eralized hypergeometric functions. 
It is easy to show [16[ that Eq(D,0) = and 



Ef F (D,0)=0 



(29) 



The Ei (and thus E c ) value for 2-sphcrium was re- 
cently reported by Seidl 15|. However, simple expres- 
sions for any D can be obtained from Eq. (|2"8l . For ex- 
ample, 



E c (2,0) =41n2-3 

4 368 

3 ~ 277t2 
64, „ 229 



E c {3,0) 
£ c (4,0) 



£ c (5,0) 



In 2 
75 375 

24 2 650112 



35 385 875tt 2 
Eq. ([28]) also yields the large- D expansion 



E c (D,0) 



1 r2 21 r3 21 r4 

-<r 5* -\ S 4 + 

8 128 512 



(30a) 
(30b) 

(30c) 

(30d) 

(31) 
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IV. HOOKIUM 



The one-electron Hamiltonian in £-hookium is 

r d 2 D — 1 d 

dr 2 r dr 



ii 



r 



and the zeroth-order wave functions are 

D 



(32) 



(33) 



fc=i 



where x^fc is the kth cartesian coordinate of electron i, 
and dfc and b k are non-negative integers. The orbitals are 
the onc-dimcnsional harmonic oscillator wave functions 



ip a {x) = V2 a alir 1 / 2 H a {x) exp(-.T 2 /2) 



(34) 



where G is Catalan's constant [18| ■ Similar remarks per- 
tain to the HF energies with odd D, such as 



4 

3 ~ 

8 


4 r 

- 1 

7T - 

8 


27 


~ 27^ 



1 + ln(8 - 4^3) 



(41a) 



3v / 3 + 61n(8-4v / 3) 



(41b) 



Eqs (|3"5)) and (j3"9"j) also yield the large- £ expansions 

3 



E 2 (D,2) 
E? F (D,2) 
E C {D,2) 



-S 2 -—d 3 
4 32 



64 



<5 2 



256 ' 
47 

256 ' 



1024 
J7_ 

1024 



(42) 
(43) 
(44) 



where H a is the ath Hermite polynomial |18| . The energy 
difference between the eigenstates are given by 



D 



e = ^2(a k +b k ) = 2i 



(35) 



k=l 



where 2n is the excitation level, i. e. the number of nodes 
in It is not difficult to show that £q(£, 2) = D and 



E 1 (D,2) = 



V2 r(f ) 



(36) 



Both £2 and £^ F can be found by direct summation 
[lH , as in Eq. The sum includes all single and 

double excitations for E 2 , but only singles for £^ F . The 
integral (^o \ r i 2 \ ^t) vanishes unless all of the a k + b k 
are even and, in that case, it is given by 



(*o\r 12 \*i) 



i r(^)r(n+j)^ 



k=l 



In this way, one eventually finds 

. 2 



E 2 (D,2) = - 

ri 



■ D-l- 
. 2 



E^(D,2) = - 



4r(f) n 

D-l \ 2 oo 
2 



E 



4r") n! n 

. 2 / n 



E 



(1/4)' 



n! n 



Aa k -b k 

r 

V7ro fe !6 fc ! 

(37) 



(38) 
(39) 



which reduce to generalized hypergeometric functions. 

£2(3, 2) has been derived by several groups [IJ Hi 
[30j but the energies for other £ have not been reported 
before. All can be found in closed form and the first few 
are 



£2(2,2) 


= 2G-7rln2 


(40a) 


£2(3,2) 


= 1 - 


- -(l + ln2) 

7T 


(40b) 


£2(4,2) 


1 

~ 4 


2G - TT In 2 + 1 - - 
4. 


(40c) 


£2(5,2) 


5 

~ 9 


-^(4 + 31n2) 


(40d) 



V. RESULTS AND DISCUSSION 

Numerical values of E 2 , £^ F and E c , for D = 2, . . . , 8 
and m = — 1, and 2 are reported in Table |U The £2 
values for helium were found by the Hylleraas technique 
described in Section II. Other results were obtained from 
Eqs (Till), O, (Ml, © and ©. 

As m increases (for constant £), although the exact 
and HF energies decrease in magnitude from helium to 
spherium and then increase from spherium to hookium, 
the correlation energies always increase. However, the 
smallness of that increase is striking; £ c is almost inde- 
pendent of to, especially for large D. The correlation 
energies of helium and hookium differ by only 8% for 
D = 2, and this drops to just 2% for D = 8. 

As D increases (for constant to), all of the energies 
decrease dramatically and the correlation energies fall by 
a fe ahT|^st[-tt^ orders of magnitude between D = 2 and D = 
8. H crrick) and Stillinger have explained this in £-helium 
[19t f by observing that the Jacobian (ri^) 15-1 creates 
a "dimensionality barrier" that keeps both electrons far 
from the nucleus and therefore allows them to avoid each 
other more easily when D is large. Similar arguments 
apply to £-spherium and £-hookium and, presumably, 
in general. 

The observed dependence of the correlation energy on 
D and to is consistent with the large- 1? expansions (|2ip . 
(13H and (l44"l). all of which take the form 



E c (D,m) ~ -<5 2 /8- CS 3 



(45) 



where the coefficient C ps 1/6 varies slowly with to. Such 
an expression implies that £ c depends primarily on the 
dimensionality of space in which the electrons move but 
with a small correction from the shape of the confining 
external potential. 

We conjecture that Eq. (|4"S")) is true for all confining 
external potentials V(r). To explore this, it would be 
useful to extend Table U to include "airium" \V(r) = 
r], "ballium" [3l| [V(r) = r°°] and other such systems. 
These studies will be reported elsewhere. 
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TABLE I. Second-order energies and limiting correlation energies in two-electron systems. 



System 



D = 2 



D = 3 



D = 4 



£> = 5 



D = i 



D = 7 



D = 8 



Second-order exact energies, —E2(D,m), from ([3} 

Helium -1 0.632740 0.157666 0.070044 0.039395 0.025208 0.017501 0.012854 

Spherium 0.227411 0.047637 0.019181 0.010139 0.006220 0.004189 0.003007 

Hookium +2 0.345655 0.077891 0.032763 0.017821 0.011153 0.007622 0.005533 

Second-order HF energies, -Ef F {D,m), from © 

Helium -1 0.412607 0.111003 0.051111 0.029338 0.019020 0.013325 0.009852 
Spherium 00000000 

Hookium +2 0.106014 0.028188 0.012904 0.007382 0.004776 0.003342 0.002469 

Limiting correlation energies —E c (D,m), from ([5]) 

Helium -1 0.220133 0.046663 0.018933 0.010057 0.006188 0.004176 0.003002 

Spherium 0.227411 0.047637 0.019181 0.010139 0.006220 0.004189 0.003007 

Hookium +2 0.239641 0.049703 0.019860 0.010439 0.006376 0.004280 0.003065 
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